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ABSTRACT 

We have made a comprehensive transit search for exoplanets down to about 2 Earth radii in the 
HD 209458 system, based on nearly uninterrupted broadband optical photometry obtained with the 
MOST (Microvariability and Oscillations of Stars) satellite, spanning 14 days in 2004 and 44 days 
in 2005. We have searched these data for limb-darkened transits at periods other than that of the 
known giant planet, from about 0.5 days to 2 weeks. Monte Carlo statistical tests of the data with 
synthetic transits inserted allow us to rule out additional close-in exoplanets with sizes ranging from 
about 0.20 - 0.36 Rj (Jupiter radii), or 2.2 - 4.0 (Earth radii) on orbits whose planes are near that 
of HD 209458b. These null results constrain theories that invoke lower-mass planets in orbits similar 
to HD 209458b to explain its anomalously large radius, and those that predict "hot Earths" due to 
the inward migration of HD 209458b. 

Subject headings: extrasolar planets: HD 209458 - planetary systems - methods: data analysis 



1. INTRODUCTION 

Since the discovery of the first exoplanet aroun d 
a Sun-like star, 51 Pegasi b ijMavor fc Quelozl fl995h . 
over 200 other ex oplanets have been discovered to date 
(|Schneider 1120061 ). Of these, a significant number have 
been "hot Jupiters" - gas giant planets orbiting ex- 
tremely near their parent stars. With the launch of 
the MOS T (Microvariability & Oscillations of ST ars) 
satellite ()Walker et all 120031 iMatthews et all 120041) in 
2003 , and the upcoming launches of t he COROT (iBaglin I 
2001 iBarge et al | [200l and Kepler (jBorucki et al.ll2004 
Basri et al. 2005) space missions, the age of transit 
searches with continuous space-based photometry is now 
upon us. These new observatories signal that the discov- 
ery of an Earth-sized planet orbiting a Sun-like star is a 
prospect that can now be seriously anticipated. 

1 Based on data from the MOST satellite, a Canadian Space 
Agency mission, jointly operated by Dynacon Inc., the University 
of Toronto Institute of Aerospace Studies and the University of 
British Columbia, with the assistance of the University of Vienna. 

2 Dept. Physics & Astronomy, University of British 
Columbia, 6224 Agricultural Road, Vancouver. 
1Z1, Canada; croll@astro.utoronto.ca, matthewsC 
rowe@phas.ubc.ca, kuschnig@phas.ubc.ca, gladmanC 



BC V6T 
phas.ubc.ca, 
phas.ubc.ca, 



ccameron@phas.ubc.ca 

3 Sumus Technology Ltd.; arwalker@sumusltd.com 

4 Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138, USA; sasselov@cfa.harvard.edu 

5 1234 Hewlett Place, Victoria, BC V8S 4P7, Canada; gor- 
donwa@uvic.ca 

6 University of California Observatories, Lick Observatory, Uni- 
versity of California, Santa Cruz, CA 95064, USA; lin@ucolick.org 

7 Department of Astronomy and Physics, St. Mary's University, 
Halifax, NS B3H 3C3, Canada; guenther@ap.stmarys.ca 

8 Obs du mont Megantic, Dept de physique, Univ de Montreal, 
CP. 6128, Succ. Centre-Ville, Montreal, QC H3C 3J7, Canada; 
moff at @ astro . umontreal . ca 

9 Dept. Astronomy & Astrophysics, David Dunlap Obs., Univ. 
Toronto, P.O. Box 360, Richmond Hill, ON L4C 4Y6, Canada; 
rucinski@astro.utoronto.ca 

10 Institut 



The first transits of an exoplanet were observed in the 
system HD 209458, whose planet was discovered from 
radial velocity measurement s by Mazeh et al. ( 2000) . 
The transits were detect ed bv lCharbonneau et alj (2000) 
and I Henry et al.l ([2000) . The discovery of this transit- 
ing planet has led to concentrated studies of the sys- 
tem and hence improved knowledge of the star's char- 
acter i stics (e.g..rW ittcnmvc r et al.ll2005l : iMandel fc Ago] 
I2002t iKnutson et all 120071) . which is critical for the de- 
termination of the radius and other properties of the 
known exoplanet and any additional transiting exoplan- 
ets found orb iting HD 209458. HP 209458 is a GO V, 
V=7.65 star (|Laughlin et alj |2005[ ) with effective tem- 
perature T e g = 6000 ± 50 K, luminosity L = 1.61L Q 
(jMazeh et all [2000) . mas s M = 1.10 ± 007Afo , and ra- 
dius R = 1.13 ± O.O2i? (jKnutson et alJl2007f >. 

There are theoretical reasons to expect smaller exo- 
planets in close orbits around HD 20945 8. The anoma- 
lously large size of the known exoplanet (|Laughlin et al.l 
2005) has led to suggestions that the planet may 
be "puffed up" due to tidal heating involving the 
interaction of a second planet (lLaughlin et al.1 120051 : 
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I Raymond et all (|2006l ) predict that the inward migration 
of hot Jupiter planets will lead to Earth-mass and super- 
Earth-mass planets in close orbits to their parent stars 
in systems such as HD 209458. 

Detection of exoplanets around other stars via 
transit measurements is a field of increasing impor- 
tance. Groun dbased wide-field photome tric surveys 
such as TrES ([Alonso et aT1l2Q0j), HAT (jBakos et all 
l2002lL OGLE (lUdalski et all 120021). and SuperWASP 
(|Street et al. II2004D are simultaneously monitoring thou- 
sands of nearby bright star s for evidence of tr ansits. 
The discoveries of TrES-1 (I Alonso et al l I2004D. XO- 
lb dMcCullough et all 120061). TrES-2 (10 'Donovan et all 
I2006D. and HAT-P-lb (TBakos et all I2006D . have illus- 
trated the feasibility of discovering extrasolar planets via 
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the transit method. Other groundbased surveys have fo- 
cused primarily on globular and o pen clusters and have 
returned a series o f null results (IWeldrake et alj 120051: 
Hidas et al l 120051: Hood et all [20051: Ivon Braun et all 



20051: iBurke et aLll2006l : iMochejska et al.ll2006l ). Space- 
based HST results returned a null-result in the globular 
cluster 47 Tucanae (Gilhlan d~eraT1l2000[) . while more re- 
cently ISahu et "all (|2006l ) detected 16 transiting planet 
candidates near the Galactic bulge. 

All groundbased transit searches are inherently lim- 
ited by time sampling: the day-night cycle, weather and 
sometimes observing time allocation. In addition, their 
photometric precision means they are sensitive to only 
planets similar in scale to our giant planets. Space-based 
transit searches from satellites in inclined equatorial or- 
bits (such as HST) also have limited time sampling. How- 
ever, space-based transit surveys from platforms with 
large continuous viewing zones (CVZ) on the sky offer 
the potent combination of excellent extended-time cov- 
erage and high photometric precision. 

In this paper, we describe such a search for transits 
in the HD 209458 system with the MOST satellite, with 
nearly continuous time coverage over two epochs span- 
ning two and six weeks in 2004 and 2005, respectively. 
The methods refined and developed here to take advan- 
tage of this precise, high duty-cycle photometry allow 
us to reliably detect transits as shallow as ~0.5 mmag, 
corresponding to an exoplanet as small as about 0.2 Rj, 
or 2 R^. These methods are also applicable to other 
space-based transit searches, such as future searches with 
MOST, as well as those planned with the COROT and 
Kepler space missions. 

In S}2] the MOST photometry of HD 209458 is de- 
scribed. The transit search technique is discussed in Ej3j 
The Monte Carlo statistics used to estimate the sensitiv- 
ity of the transit search are specified in The transit 
search routine is applied to the MOST HD 209458 data 
set in [JSJ and the results are presented in ^ including 
a discussion of the impact of these results on theories 
relating to other putative exoplanets in the system, and 
the characteristics of the known exoplanet HD 209458b. 

2. MOST PHOTOMETRY OF HD 209458 

The MOST satellite was launc hed on 30 June 2003 and 
its initial mission is described bv lWalker et~aT] (|2003t ) and 
iMatthews et~all (l200l ). A 15/17.3 -cm Rumak-Maksutov 
telescope feeds two CCDs, one originally dedicated to 
tracking and the other to science, through a single, cus- 
tom, broadband filter (350 - 700 nm). MOST was placed 
in an 820-km circular Sun-synchronous polar orbit with 
a period of 101.413 minutes. From this vantage point, 
MOST can monitor stars in a continuous viewing zone 
(covering a declination range +36° < S < —18°) within 
which stars can be monitored without interruption for 
up to 8 weeks. Photometry of very bright stars (V < 6) 
is obtained in Fabry Imaging mode, in which a Fabry 
microlens projects an extended image of the telescope 
pupil illuminated by the target starlight to achieve the 
highest precision (see Matthews et al. 2004). Fainter 
stars (down to about V ~ 12) can be observed in Direct 
Imaging mode, in which defocused images of stars are 
monit ored in Science CCD subrasters (see iRowe et al.l 
I2006al) . 



MOST observed HD 209458 in Direct Imaging mode 
over two sep arate epochs: for 1 4 days during 14 - 30 
August 2004 (|Rowe et al.ll2006al). and 44 day s during 2 
August - 15 September 2005 (jRowe et al.ll2007[ ). The first 
run was a trial, and testing of new onboard software led 
to several data interruptions, reducing the overall duty 
cycle of the raw photometry to about 85%. The sec- 
ond run returned data without significant interruptions, 
achieving an overall duty cycle of the raw photometry 
of 97%. The exposure time of the observations was 1.52 
sec, and the sampling rate was 10.0 sec for both epochs 
of data. Approximately 112,000 and 361,000 individual 
observations were taken in the 2004 and 2005 epochs, 
respectively. 

The reduction of the raw photometry downloaded from 
the satellite (and converted into FITS files) for both 
epochs was performed by JFR. The reduction is similar 
to that applied to groundbased CCD photometry, but 
is nondifferential, and incorporates both aperture and 
PSF (point spread function) fitting. It corrects for cos- 
mic ray hits (especially frequent during satellite passages 
through the South Atlantic Anomaly [SAA]), the vary- 
ing background due to scattered earthshine modulated 
at the satellite orbital period, and flatfielding effects. 
T hose data reduction tech niques used are fully described 
in lRowe et all ()2006a| ) and lRowe et~all (|2007f i. 

2.1. Additional filtering and selection of the data 

For this transit search, additional filtering and selec- 
tion of the photometry was done to optimize the data for 
the application of Monte Carlo statistics. The filtering 
is an automated part of the transit search routine, and 
is thus briefly summarized in ij3.ll This filtering step is 
described in detail here. 

Data obtained during passages through the SAA are 
conservatively excised from the light curve, due to the 
increased photometric scatter during those MOST or- 
bital phases, without reducing seriously the phase cov- 
erage of the exoplanetary periods searched. The tran- 
sits of the known giant planet HD 209458b are removed, 
at the orbital period P — 3.52474859d determined by 
iKnutson et al.l (|2007f ). Because of the modulation of 
stray earthshine with MOSTs 101.4-min orbital period, 
the data were phased to that period, and segments show- 
ing the most noticeable effects from stray light were also 
removed. The coverage in phase per 101.4-min MOST 
orbit following this cut was 63% and 58% for the 2004 
and 2005 epochs, respectively. The stray light back- 
ground can also be modulated at a period of 1 day and 
its first harmonic (due to the Sun- synchronous nature of 
the MOST satellite orbit), so sinusoidal fits with periods 
within 1% of 1 and 0.5 d were subtracted from the data. 
After these cuts, any remaining outliers greater than 6a 
were excised. This sigma cut removed very few data- 
points, as few points were such extreme outliers. The 
magnitude of the injected transits (Q was always at a 
level much less than this sigma-cut. The data were also 
median-subtracted, for reasons outlined in §3.51 The au- 
tomatic filtering and selection step removed 45% of the 
original data (mostly due to the SAA and stray light 
corrections). The root mean square of the data following 
this additional filtering and selection step is 0.0035 mag. 

The resulting HD 209458 light curves are plotted un- 
binned, and binned in 30-min intervals in Figure [TJ 
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These filtered data (and the original reduced data) can 
be downloaded from the MOST Public Data Archive at 
www.astro.ubc.ca/MOST. 

3. TRANSIT SEARCH ALGORITHM 

Our transit search routine has been developed and re- 
fined to take specific advantage of the unique time cov- 
erage of the MOST data. For HD 209458 in particu- 
lar, it is designed to make use of the well determined 
stellar and exoplanetary orbital parameters of the sys- 
tem. The search is intended to run automatically with- 
out user intervention to facilitate the generation of Monte 
Carlo statistics of the detection thresholds and their cor- 
responding significances. 

The entire routine (described in the following subsec- 
tions) includes: ( <j3 . 1 [) an automatic filtering and selec- 
tion step; f H3.2[) a transit search algorithm, adapted from 
the EEBL S (Edg e Effec t Box-fitting Least Squares) algo- 
rithm of lKovacsl (|2003h sensitive to box-shaped transits; 
( H3.3|) selection criteria, to narrow the EEBLS candidates 
to a select few that can be tested for astrophysical plau- 
sibility; f §3.4p transit model fitting, to identify realistic 
limb-darkened transit candidates for exoplanets of vari- 
ous radii, orbital periods and inclination angles; and fi- 
nally ( §3.5|) detection criteria, to differentiate between 
likely transiting candidates and false positives. This au- 
tomated search routine can also be applied to other tran- 
sit search experiments in which the characteristics of the 
target stars are well known or inferred. 

3.1. Automated filtering and selection scheme 

The initial reduction scheme of JFR as described in 
fj2]is almost entirely independent of the magnitude of an 
individual data-point. Instead other indicators such as 
the sky background level are used to remove or adjust the 
data, or the associated uncertainty of an individual data- 
point. Therefore, inserting synthetic transits in the light 
curve after this initial reduction to produce Monte Carlo 
statistics is a valid procedure. To ensure proper appli- 
cation of those Monte Carlo statistics for this and other 
data sets, we incorporate an automated data filtering and 
selection scheme in the transit search to filter out stray 
carthshine, and to remove additional outliers and the 
dominant transits of the known exoplanet. This process 
was described in detail in H2.ll It is briefly summarized 
here to ensure that it is obvious that this process is an 
automated part of the transit search routine applied to 
the Monte Carlo statistics to come. The thresholds and 
amplitudes used in this automatic filtering scheme have 
been optimized for their application to the MOST HD 
209458 data; obviously these thresholds and amplitudes 
would have to be similarly optimized for their application 
to other MOST data-sets. 

3.2. EEBLS (Edge Effect Box- fitting Least Squares) 
algorithm 

To search for transits in the HD 209458 system at peri- 
ods other than that of the known exoplanet HD 209458b, 
we use a slight modificat ion of the Box-fitt ing Least 
Squares (BLS) algorithm (|Kovacs et al.l I2002T ). namely 
the Edge Effect Box-fi tting Least Squares (EEBLS) al- 
gorithm (jKovacs Il2003h . The EEBLS routine flags puta- 
tive transit candidates by searching for box-shaped tran- 
sits in a light curve. The signal strength of a putative 



transit is indicated by the Signal Residue (SR), as de- 
fined in lKovacs et~aTT(|2002f ). The SR is roughly equiv- 
alent to the x 2 for a box-shaped tr ansit as indicated in 
the appendix of IBurke et ail (|2006fl . The EEBLS rou- 
tine returns the period, phase, and drop in magnitude 
of the supposed transit. The EEBLS algorithm was 
chosen becaus e it has been used extensively in similar 
searches (e.g. IBurke et al.l 120061: iMocheiska et al.l 120061: 
ISahu et al.ll2006h . and quantitative com parisons indica te 
it is as good as others in the literature (|Tinglevll2003l ). 

As is the case with all computational routines that 
must handle a large number of data and an extensive 
search parameter grid, sufficient resolution must be bal- 
anced against computational duration. To decrease the 
latter, we employ even logari thmic period spacin g of the 
search grid as suggested by IBurke et al.l (|2006l ) . With 
even logarithmic period spacing, the subsequent period 
searched, P2, is related to the previous, Pi, by P2 = 
Pi(l + 77), where r\ <C 1. Although even frequency 
spacing was used in the original formulation of BLS by 
iKovacs et all (2002). even logarithmic period spacing re- 
tains high sensitivity to transits with both short and long 
orbital periods without increasing computational dura- 
tion. Both even frequency and even logarithmic period 
spacing are valid with sufficient resolution. We chose the 
value 7] = 0.000168 to attain the especially fine period 
resolution needed to resolve putative transits occuring 
in the two epochs of MOST data separated by the near 
one-year gap. 

In this work, we search the data set for periods greater 
than P m i n = 0.5<i, and less than P max — 14. 7d. The 
maximum period limit is set by our requirement to ob- 
serve at least 3 hypothetical transits of a putative planet 
in the longer 2005 run (44 days). The minimum period 
was set to be 0.5d, likely below the astrophysically rea- 
sonable period of the orbit of a stable exoplanet in the 
HD 209458 system, as it would correspond to a semi- 
major axis less then ~2.5 times the stellar radius. 

The EEBLS routine was also modified to enhance its 
detection sensitivity. For each trial period, the EEBLS 
routine was changed to record the most significant box- 
car dimming (transit) and box-car brightening signal in- 
dependently. The original formulation recorded the most 
significant box-car event, regardless of whether it was a 
dimming or a brightening. This resulted in a modest de- 
crease in sensitivity as the noise floor increased because 
the algorithm was automatically including brightening 
events that could not possibly be transits. This reformu- 
lation allows one to better differentiate possible transits, 
from systematic brightening events. 

In addition the minimum and maximum fractional 
transit lengths, Qmi and Qma respectively, that the 
EEBLS algorithm searches each trial period for were 
set to be variable, rather than constant as in the orig- 
inal formulation. Th e small planet approximation of 
iMandel fc Ago] (|2002D , as discussed below, was used to 
determine the approximate maximum fractional transit 
length for each trial period, Qmp. For this approxima- 
tion a planet with radius R p = 0.4Pj, and with inclina- 
tion angle i = 90° was used to determine this maximum 
fractional transit length. This change is significant be- 
cause the difference in the maximum fractional transit 
length falls from nearly Qmp w 0.13 at the minimum 
period of 0.5d, to Qmp « 0.015 at the maximum period 
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TABLE 1 

EEBLS INPUT PARAMETERS 



vax 


definition 


value 


Np 


Number of Period Points searched 


20 000 


V 


logarithmic period step 


0.000168 




min. period submitted to EEBLS 


0.5d 


Pmax 


max. period submitted to EEBLS 


14.7<2 


Qmi 


min. fractional transit length 


0.75 xQmp 


Qma 


max. fractional transit length 


l.lxQmp 


Nb 


Number of phase bins 


20.0/Qmi 



of 14. Id. This change allows one to thus use the known 
characteristics of the star that is investigated to precisely 
tailor the EEBLS search routine. Qmi and Qma were set 
significantly lower, and slightly larger, respectively, than 
this maximum fractional transit length for each period; 
that is Qmi was set as 0.75 times as large, and Qma was 
set as 1 . 1 times as large as the maximum fractional tran- 
sit length for each period. These values of 0.75 and 1.1 
times as large as the maximum fractional transit length 
were chosen because they span the parameter range of 
interest from a large edge-on planet to a small grazing 
planet. 

The final change that was implemented was to set a 
variable number of bins in the folded time series at each 
test period, Nb, rather than the constant number as in 
the original formulation. This was again related to the 
fact that the maximum fractional transit length varies 
significantly from the minimum to maximum period in- 
vestigated in this formulation. The decision was made 
to ensure that there was a constant number of bins in 
the folded time series during the in-transit component 
for the smallest fractional transit length, Qmi, for each 
trial period. This constant number of bins was set ar- 
bitrarily at 20, as it was believed this was large enough 
to properly determine the phase of a putative transit- 
ing planet. That is for each trial period Nb was set 
as 20.0/Qmi. Nb thus varied from as low as «200 
for the minimum period of 0.5d, to as high as ~1570 
for the maximum period of 14. Id. These changes not 
only increase the sensitivity of the algorithm to long- 
period events, but also marginally decrease the compu- 
tational duration of the EEBLS algorithm. The EE- 
BLS input parameters are summarized in Table [1] The 
modified version of this algorithm in C is available at: 
www.astro.ubc.ca/MOST/EEBLSmodified.cpp 

3.3. Transit Selection Criteria 

The EEBLS algorithm returns the Signal Residue (SR) 
versus period for each of the Np periods investigated. 
Selection criteria have been developed to pick the most 
likely genuine transit candidates. These selection criteria 
are described here. Of all the periods, only those that 
are strong local maxima in SR versus period are flagged 
for further investigation. A strong local maximum is de- 
fined as one for which SR is greater than the 0.003 x 
nP neighbouring period points in both directions. Of 
the strong local maxima, the candidates with the high- 
est 2 SR signals are chosen for further analysis. As we 
often noted that simply picking the top SR values bi- 
ases one's selection against short-period planets of small 
radius, a second selection criterion was implemented to 



ensure adequate sensitivity to short-period events (~0.5 
< P <~4.0). The top 2 candidates with SR values 1.8x 
greater than the SR noise floor were also selected for ad- 
ditional study. The SR noise floor is defined as the mean 
of the SR values for all period points within 0.2 x nP of 
the period of interest. These two selection criteria nar- 
rowed the list of prospective candidates to a maximum of 
four of the most likely transit candidates. Extensive sta- 
tistical tests (described in §4]) both verified and refined 
the above criteria. 

The EEBLS SR spectrum for the MOST HD 209458 
2004 and 2005 light curves is shown in Figure [2j the can- 
didates that passed the selection criteria for these light 
curves are also identified in Figure^ 

We also optimize the determination of the exact period, 
P, and phase, (ft, for these candidates by rerunning the 
EEBLS algorithm with greater period, and then phase 
resolution, in the immediate period space of the candi- 
date event. First, we run the EEBLS algorithm with 
greater period resolution from a minimum to maximum 
period of 4 period points less to 4 period points greater, 
respectively, than the best-fit period. An order of mag- 
nitude more period points are then used (100 points), 
to more accurately determine the period, at which the 
greatest SR signal is returned. The period that returns 
the highest EEBLS SR is then investigated with the EE- 
BLS algorithm at this one specific period with 5 times as 
many phase points as previously. These two steps are of 
negligible computational duration, but considerably in- 
crease the accuracy with which the period, and phase are 
returned. 

3.4. Realistic Transit Fitting 

Since the discovery of the giant exoplanet HD 209458b, 
there has been intensive study of the star HD 209458, 
such that its basic characteristics are now w ell known 
(jWittenrnver et all 120051 : iKnutson et aT1l2007f ). This al- 
lows for considerable improvement in the shape, dura- 
tion, and depth of the transit model as opposed to the 
original box-shaped transit used in the EEBLS routine. 
Specifically, we search for realistic limb-darkened transits 
of a specific period, P, planet radius, R p , orbital inclina- 
tion angle, i, and phase of tran sit, (ft. Our model uses the 
small planet approximation of iMandel &: Agol (|2002f ) to 
set the shape, length and depth of the transit. T he small 
planet approximation of IMandel &: Ago] (|2002j ) is valid 
for -g^ < 0.1, and thus should be valid for HD 209458 for 
all transiting planets below Rp « 1.1 Rj. We have only 
searched for transiting planets whose orbits are circular 
(e = 0).; The impact of eccentricity is discussed in 52J 
The HD 209458 system parameters used for our analysis 
are given in Tabled 

The transit candidates are then fitted by our limb- 
darkened transit model. For the period P returned by 
the EEBLS algorithm, two inclination angles i are inves- 
tigated. One inclination angle is edge-on (i = 90°), while 
the other is evenly spaced between the minimum inclina- 
tion that would produce a transit (cosi ~ R*/a, where a 
is the semi-major axis), and an edge-on (i = 90°) orbit. 
For short period events, where the minimum inclination 
that would produce a transit is below 78°, the minimum 
inclination for the initial guess is arbitrarily set to 78°, 
as it was found this increased the detection efficiency. 
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TABLE 2 

System parameters 



var definition value 







Stellar radius 


1.125 Rq a 


M, 




Stellar Mass 


1.101 Mq a 


Cl 


non- 


-linear limb-darkening parameter 1 


0.410769 b 


C2 


non- 


-linear limb-darkening parameter 2 


-0.108909 b 


C3 


non- 


-linear limb-darkening parameter 3 


0.904020 b 


C4 


non- 


-linear limb-darkening parameter 4 


-0.437364 b 



a Parameters obtained from IKnutson et al.l |(2003) b These 

limb-darkening para meters are fully described in 

Mandcl & Agol (2002) and obtained from Knutson ct al. 
H2007I) and lRowe et all lt20O7T) 



For example a period of approximately 4 days will have 
a minimum inclination angle of 84°, and thus the two 
guesses will be 87°, and 90°, respectively. For these incli- 
nation angles, the planet radius R p , phase 0, and inclina- 
tion i, of a putative transit are fit through a Marquardt- 
Levenberg non-linear least squares algor ithm (jjVlarquardt 
U963HLevenberdll944t|Press et al.|[l99l) . Different start- 
ing values for the orbital inclination are used, since it 
was found through Monte Carlo analysis of simulated 
transits (Q that this improved the efficiency in return- 
ing the correct value of i. The initial guesses in radius, 
R p , and phase, 0, given to the Marquardt-Levenberg al- 
gorithm are the characteristics returned by the EEBLS 
routine, where the dimming in brightness is translated 
into the size of a hypothetical edge-on planet using the 
known stellar radius. Limits were placed on the range 
of parameters that the Marquardt-Levenberg algorithm 
could explore, away from the starting points returned by 
the EEBLS algorithm, to ensure that the fit would re- 
main in the parameter space that was flagged as signifi- 
cant by the EEBLS algorithm. These were AR p w 0.10, 
arccos(i?*/a) < i < 90.0°, and A</> « 0.05. Fitting for 
the period as well using the Marquardt-Levenberg non- 
linear least squares algorithm was not attempted here, 
as the period optimization discussed above was deemed 
sufficient. The values of R p , i and (f> that produce the 
minimum \ 2 value are recorded as the best transit can- 
didate for that period. 

An attempt was made to bypass the EEBLS algorithm 
altogether and solely use the above fitting method to 
search for realistic limb-darkened transits in the MOST 
HD 209458 data. However, this approach proved very 
computationally intensive, and so we adopted and recom- 
mend the compromise described above, which achieves 
sufficient accuracy without significant increase in com- 
putational duration. 

3.5. Transit Detection Criteria 

We have developed quantitative criteria to differentiate 
bona-fide transiting planet candidates from false positive 
detections, without the need for detailed visual inspec- 
tions of every set of candidate events in the light curve. 
We compare each modeled transit candidate to an anti- 
transit (brightening) model, and a constant brightness 
model. The latter assumes a star of constant magnitude, 
given by the mean of the light curve. The transit model 
consists of an out-of-transit component and an in-transit 
component (modeled as described in 33. 4p . The out-of- 



transit component has a constant brightness level given 
by the median of the light curve; the median was cho- 
sen, as opposed to the mean, because it should remain 
relatively unchanged regardless of whether one or more 
other transiting planets are present in the HD 209458 
system. During transit a realistic limb-darkened transit 
of a spherical planet is used with parameters R p , i, and 
4>, as returned by the fit described in 33.41 and the small 
planet approximation of iMandel fc Agol (|2002l ). 

The statistic that has been used to quantify the im- 
provement of the transit model over the constant flux 
model is Ax 2 %, the percentage improvement in \ 2 of the 
transit model (Xt) over the constant flux model (xc)- 
That is A X 2 % = 100 x x o~xl _ The ant i-transit model 
is identical to the transit model except that the transit 
causes brightening as opposed to dimming. The best-fit 
anti-transit model was determined through the identical 
method as outlined in steps §3.21 through §3.41 Simulta- 
neously determining the best-fit transit, and anti-transit 
model results in a minor increase in computational du- 
ration, but this is a justified sacrifice as it provides an 
independent method for determining the detection sensi- 
tivity. In a similar fashion, the improvement of the anti- 
transit model {x'at) over the constant-brightness model 
is given by Axt% = 100 * x c~*^ . The quantitative 

detection criteria were motivated bv lBurke et~afl (2006) 
and have been refined through the Monte Carlo tests de- 
scribed in §3 

A believable transit was defined as one that satisfied 
the following criteria: 

• Of the transits selected by the selection criteria 
( §3.3|) the transit of interest showed the greatest 
improvement in A% 2 %. 

• The ratio between the improvement of the transit 
model over the constant-brightness model versus 
the anti-transit over the constant-brightness model 
is at least 1.01 (A x 2 %/A x l% > 1.01). 

• The quantity /, as defined bv iBurke et alj (|2006f ) 
and explained below, is less than 0.65 (/ < 0.65). 

The first criterion is intended to be conservative, as it 
helps to avoid identifying harmonics of the period of an 
obvious tr ansit candidate. Th e second criterion was sug- 
gested by IBurke et al.l (2006) for a groundbased transit 
survey, and is primarily used as a threshold to differenti- 
ate astrophysically credible transit candidates from false 
positives. We use a much more sensitive value for the 
ratio (Ax 2 %/Axt% > 1-01), due to the increased sen- 
sitivity of the MOST space-based photometry, as well as 
the results of the Monte Carlo tests of §3 The level that 
has been chosen requires the best-fit transit to be only 
slightly more prominent than the best-fit anti-transit 
candidate. It is important to note that the best-fit anti- 
transit is not restricted to at the same period as the 
best-fit transit. The best-fit anti-transit is determined 
in a completely analogous fashion to the best-fit tran- 
sit, and thus the parameters of the best-fit transit and 
anti-transit are often completely independent. The moti- 
vation for selecting this sensitive value of 1.01 was based 
on the experience gained from the Monte Carlo results, 
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as it was the set at that level at which a transit could 
be reliably recovered from the data-set without inducing 
significant numbers of false-positives. The full explana- 
tion of how the value of Ax 2 %/Ax?_% > 1.01 was de- 
riv ed is given in flU T he rationale for this cut, as noted 
bv iBurke et al.l (|2006h . is that in the case that the light 
curve shows periodic variability or significant brightening 
or dimming trends, a transit search will often assign con- 
siderable significance to these events, possibly resulting 
in false positive detections. Although this cut has pri- 
marily been designed for wide-field, or open and globular 
cluster transit searches, where eliminating targets that 
show obvious sinusoidal variation is desirable, it is still 
applicable to the present search to account for possible 
variations and trends due to stray earthshine or intrin- 
sic variability in the star HD 209458. It should also be 
noted that only improvements in the anti-transit model, 
A%?_% (Ax-% > 0), are accepted. This is because, in 
the limit of extremely deep transits, on the order of those 
visible by eye such as those produced by HD 209458b, the 
anti-transit statistics are forced to return worse values of 
Ax-%, as no significant brightening event is detected. 
This is because the median rather than the mean is used 
for the out-of-transit component of our transit model. 
This effect was not observed in the current application, 
or in the Monte Carlo tests, but it is a useful caveat for 
other transit searches. 

The third criter ion invokes the parameter /, which 
IBurke et al.| (|2006Q defines as f=Xkth/Xtotal- Here xl 
refers not to the typical x 2 , but rather to summing the 
following quantity for the k th transit: xtth = ( m i/ a i)i 
where and <7j are the magnitude and uncertainty, re- 
spectively, of the i th measurement during the k th transit. 
The criterion requiring / < 0.65 roughly corresponds to 
observing the transit at least one and a half times, assum- 
ing similar noise. This criterion is useful in eliminating 
long-period false-positive events, as it ensures that the 
dimming behaviour that caused the EEBLS routine to 
flag the event as highly significant, is actually observed at 
least one and a half times. This cut is of only limited sig- 
nificance to our HD 209458 search, given the sensitivity, 
time span, and high duty cycle of the MOST photome- 
try. However it remains a safeguard for some long-period 
events in our current search and for future ground and 
space-based searches. 

Avoiding false positives is necessary so as to not place 
unjustifiably sensitive limits in the Monte Carlo statis- 
tics of SI These or similar criteria should be suitable 
for other ground and space-based transiting exoplanct 
searches. 

4. MONTE CARLO STATISTICS 

To assess the sensitivity of the aforementioned search 
routine and the MOST data set to other transiting plan- 
ets in the HD 209458 system, simulated transits for plan- 
ets of various radii and orbital parameters were inserted 
into the MOST photometry and Monte Carlo statistics 
of the transit recovery rate were generated. 

Realistic limb-darkened transits due to planets with 
various radii R p , orbital phases <fi, periods P, and 
inclinations i, were inserted into the 2004 and 2005 
MOST HP 209458 data . The small planet approxima- 
tion of iMandel fc Ago! (|2002t) was used. The differ- 
ences between the more accurate non-linear model of 



IMandel fc Ago] (|2002f ) and the small planet model ap- 
proximation are negligible for the size of planets inserted 
in this data-set, even for transits with short periods, and 
thus semi-major axes only a few times the stellar ra- 
dius. These modified data were then subjected to the 
analysis described above. Transits were inserted with 
logarithmic period spacing (as discussed in ij3.2p . with 
Vinp — 0.095 in the period range 0.55<i < Pi„ p < 14. 5d 
In total, 37 period steps were used for the 90 and 88° 
inclination angle cases. For other inclination angles tran- 
sits were inserted with logarithmic period spacing until 
the period exceeded the maximum period, or semi-major 
axis, that would produce a transit, a < (i?* + Rp)/cosi. 
For each trial period, simulated transits corresponding 
to 9 different exoplanet radii were inserted, sampling the 
period-radius space of interest. For each period and ra- 
dius, 110, 65, or 25 phases were inserted as summarized 
in Figures [4] and [6l For each of these points, the phase 
4> was generated randomly to be in the range < 4> < 
1. Because a 2.4 Ghz Pentium processor with 1 Gbyte of 
memory can perform the transit search algorithm on an 
individual MOST HD 209458 data set in -10 minutes, 
exploration of the entire grid just mentioned for all in- 
clinations involves ~ 10 5 iterations and 2.0 CPU years. 
The calculation was performed on the LeVerrier Beowulf 
cluster in the Department of Physics and Astronomy at 
the University of British Columbia using 45 dual-CPU 
compute nodes. 

An inserted transit was judged to be detected if the 
parameters ip and P returned by the transit search al- 
gorithm were sufficiently close to the input values, (f>i np , 
Pinp- The returned period had to satisfy the following 
criteria: \P - P inp \ < O.ld and \P/P mp - 1| < 1%. The 
limits of the criterion for (f> were dependent on the orbital 
period, because as the period of the putative transiting 
planet decreases, the fractional transit length increases 
accordingly. Thus the accuracy required in the determi- 
nation of <p was relaxed for shorter periods. The criterion 
on (f> is: \<f> - (j) mp \ < 0.09 - 0.0054(P inp - 0.5d)/d. Obvi- 
ous multiples of the period of the inserted planet, up to 
4 times the inserted period, Pi np , as well as half-period 
(P m \Pinp), and one-third period (P « \Pinp) solu- 
tions as returned by the transit search algorithm (fj3) 
were also accepted, as it was found for a low percent- 
age of cases that harmonics or subharmonics of the in- 
serted trial period were flagged as the best candidates. 
For the 90-degree inclination angle the best candidates 
had a period double, triple, or quadruple the inserted pe- 
riod approximately 1.0%, 0.5%, and 0.3%, respectively, 
of the total transit detections in the Monte Carlo analy- 
sis, while half and one-third period solutions accounted 
for 0.3% and 0.1% of the aforementioned total, respec- 
tively. 

These Monte Carlo results were also used to determine 
the level above which a transit recovered from the data- 
set could be considered significant. In deriving this value 
one must be careful to closely balance the desire to use 
as low, and thus sensitive, a limit as possible with the ne- 
cessity of limiting the number of false positive detections. 
Thus, to satisfy this rationale the following method has 
been used to determine this value of Ax 2 %/A^ 2 _% that 
will denote a believable transit. To determine a suitable 
value of the transit over the anti-transit model ratio, 
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Ax 2 %/Ax 2 _%, all Monte Carlo cases where the candi- 
date returned by the transit routine was not the candi- 
date that was inserted were investigated. These events 
where the candidate returned by the routine was not the 
inserted transit will henceforth be referred to as spurious 
events. 

By initially investigating the Monte Carlo statistics 
without the criterion that Ax 2 %/Ax-% > 1.01 we were 
able to investigate the maximum efficiency of the cur- 
rent transit routine. By comparing the efficiency of 
the Monte Carlo results with and without this thresh- 
old one can quantitatively investigate the impact of this 
criterion. We restrict the investigation of these spuri- 
ous signals to those that the routine has a reasonable 
chance of actually detecting. This area was determined 
to be the level above approximately the 25% contour 
of the Monte Carlo statistics when this threshold cri- 
terion was not used. Thus for the generation of spu- 
rious signal statistics we only use those Monte Carlo 
cases with an input radius above the limit given here: 
Rp inp > [0.125 + (0.1/14.7) (P in - 0.5days))/days Rj. 

A histogram of these spurious events versus 
Ax 2 %/Ax?_% for all inserted candidates above the 
radius threshold discussed can be seen in Figure [31 It 
was decided that a suitable level of Ax 2 %/Ax?_% would 
be one that ruled out 95% of these spurious events. 
Thus the threshold level of A X 2 %/Ax 2 _% > 1.01 has 
a transit over anti-transit ratio greater than 95% of 
these spurious signals. That is, if there is a transiting 
planet observable in the HD 209458 system, other than 
the known exoplanet HD 209458b, and the routine 
returns a putative transit that exceeds the threshold 
Ax 2 %/Ax^% « 1.01, one can be at least 95% certain 
that an actual transiting planet has been recovered, 
rather than a spurious signal. Since this measure has 
been calibrated to this MOST data set, the length of 
the time series has already been accounted for. 

The fractions of times that the artificially inserted 
transits were recovered from the data for various radii, 
periods, and inclination angles are given in Figured] Also 
indicated in Figure|4]is the 68% contour limit that would 
be placed without the Ax 2 %/Ax^% > 1.01 criterion. 
The close agreement between the 68% contour with or 
without this criterion (Ax 2 %/Ax?_% ~ 1.01) indicate 
that this criterion does not significantly impact the sen- 
sitivity of the Monte Carlo statistics generated, while 
providing a robust limit so as to avoid false positives and 
spurious detections. 

We also explore the distribution of spurious signals 
across period-radii space. These values are displayed in 
Figure O and indicate that spurious signals are negligi- 
ble (below 1%) other than the intermediate areas where 
the routine has a small to moderate chance to correctly 
determine the inserted transit. In these intermediate ar- 
eas in period-radius space the inserted transit is small 
enough that correctly identifying the period, radius, and 
phase of the inserted transit using this routine is not 
guaranteed. However, the inserted transit is of sufficient 
magnitude to significantly impact the shape of the light- 
intensity curve and lead to spurious transit signals rising 
above the threshold, Ax 2 %/AX-%> > 1.01, for a small, 
but notable, percentage of cases. Accurately determin- 
ing the frequency of true false positives - where there is 



not another transiting planet in the data-set - as opposed 
to spurious signals is not attempted in this application, 
as it is difficult to produce a synthetic data-set that ac- 
curately reproduces the systematics of the real data-set. 
These systematics may be the result of imperfect removal 
of stray earthshine from the MOST data, or may, in fact, 
be intrinsic to the star. 

We also investigated the performance of our routine 
at harmonics and subharmonics of the periods of the si- 
nusoids that were removed from the data (P=1.0d, and 
P=0.5d). Due to the fact that the fractional length of a 
transit is comparatively large for short periods of 0.5 - 1.0 
d it was found to be difficult to properly recover a tran- 
sit signal near the harmonics and subharmonics of these 
periods for all cases. As expected and shown in Figure [6l 
the routine performs minimally worse near these values 
as a portion of the transit signal is removed when the si- 
nusoidal fit is removed in £13.11 As can be seen in Figure 
O however, the period range demonstrating these drops 
in survey sensitivity is very narrow (A£ « 2%), and the 
effect is very minor. Therefore very sensitive limits can 
still be set at these periods. 

The harmonics and subharmonics of th e period of the 
know n planet HD 209458b, P w 3.5247 (jKnutson et all 
|2007| ). are also expected to show decreased sensitivity 
to transits due to the fact that the transits of the known 
planet were excised from the data. This reduction of sen- 
sitivity at these periods is unfortunate considering that 
it is expected that low- m ass planets show a preference 
for th ese resonant orbits (|Thommes 1 120051 iZhou et all 
2005). The fractional length in phase that was removed 
at P ss 3.5247 for the transit was 4.1% of the total. 
Thus for the subharmonics of the known planet (P ~ 
7.05, 10.6, 14. Id) in approximately 4.1% of the cases 
the transits will coincide with the transits of the known 
planet, and thus the routine will be unable to recover 
the transiting planet, because this data will have been 
completely removed. For these periods the 99% and 95% 
contours of Figure [4] would likely be seriously affected, al- 
though the 68% contour should only be slightly affected. 
For the harmonics of the orbital period of HD 209458b, 
the situation is better, since not all of the transit events 
would be removed from the data. For periods near the 
first harmonic (P w 1.76c?), every second occurrence of 
a transit will be lost in only 8.2% of the cases. For the 
second harmonic (P s=3 1.17d), every third transit would 
be missed in 12.3% of the cases. For periods near these 
values the sensitivity limits can be expected to be only 
slightly degraded from those shown in Figure [4] Even 
for exact harmonics the routine's sensitivity should be 
only marginally worse than the limits quoted in Figure 
HI as the remaining transits should allow the correct pe- 
riod (or a multiple of the period) to be recovered. These 
resonances have already been sensitively examined using 
the transi t-timing technique on the known exoplanet HP 
209458b (lAgol fc Steffenl 120071 IMiller -Ricci et alJ l2006h. 
iMiller-Ricci et al.1 (|2006[ ) have applied this technique to 
the MOST photometry of HD 209458 and are able to rule 
out sub-Earth-mass planets near the inner harmonics of 
HD 209458b, and Earth and super-Earth mass planets 
near the outer subharmonics out to the 4:1 resonance. 

The detection limits presented here are valid for cir- 
cular orbits, but are largely applicable to orbits of other 
eccentricities. Eccentric orbits will increase, decrease, or 
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have little effect on the fractional length of the transit, 
depending on whether the transit occurs near apastron, 
periastron, or in between, respectively. Obviously, the 
cases where the fractional transit length is increased will 
favourably affect the detection limits presented, while the 
opposite cases will adversely affect the detection limits. 
For low-eccentricity orbits this effect should be negligi- 
ble, and thus the detection limits presented here should 
be largely applicable to these orbits. For non edge-on 
inclination angles, highly eccentric orbits will adversely 
affect the limits presented, as a greater fraction of the 
time the planet will not transit the star along our line 
of sight. For edge-on, high-eccentricity orbits the lim- 
its will likely be modestly, adversely affected due to the 
cases where the transit occurs near periastron. Numeri- 
cal simulations indicate that we would expect hot Earths 
or other exoplanets to have negligi ble eccentricities in th e 
period range we are investigating (|Ravmond et al.ll2006[) . 

Thus the search routine described above in should 
be able to detect planets with radii greater than the lim- 
its given in Figure 0J In the most optimistic case of an 
edge-on 90° inclination angle transit, for planets with 
periods of approximately a half-a-day to two weeks this 
limit is approximately 0.20 to 0.36 Rj, or 2.2 to 4.0 i?®, 
respectively, with 95% confidence. If one assumes a mean 
density of p ~ 3000 kg m~ 3 - a value averaging the vari- 
ous extrasolar super-Earth models discussed below - their 
respective masses would be 0.019 and 0.011 Mj, or 6 
and 35 M®, respectively. If these putative planet were 
gaseous, and thus possible hot Neptune analogues (p « 
1600 kg m~ 3 ), a planet of Rp = 4.0 i?® would result 
in a planetary mass of 19 M®. Thus the mass-period 
parameter space we have ruled out in this study, assum- 
ing these hypothetical mean densities, sets even tighter 
limits than those set by radial y elocity observations of 
the system. lLaughlin et al.l (|2005[ ) for instance ruled out 
other planets in the HD 209458 system with Msim > 
0.3M,/ for P < lOOd The planets we are thus able to 
rule out are one-fifteenth to one-third the mass of these 
radial velocity limits for edge-on transits. Full discussion 
of the mass of planets that have been ruled out in this 
study is given in H6.1[ 

5. HD 209458 TRANSIT SEARCH 

The MOST HD 209458 2004 and 2005 data-sets were 
submitted to the analysis outlined in Sj3] There were 
no transiting planet candidates that met the detection 
criteria as outlined in i )3.5l The details of the candidate 
with the greatest improvement in y 2 are given in Tabled 
and are shown in Figure [7J As this transiting candidate 
did not meet the detection criterion, (Ax 2 %/Ax'l% w 
0.893 < 1.01), as outlined in ^3.5l we do not report it as a 
putative transiting candidate. This event, with a period 
of, P k, 14.305(i, has modest statistical significance, and 
thus we report it is a possible, but unlikely, candidate. 
The characteristics of this planet are P = 14.305c!, R p 
= 0.261i?j (2.92i?®), and i = 88.8°. If this planet is a 
super-Earth (p w 3000 kg m -3 ) the mass of this putative 
planet would be approximately 13.6 Af® (0.043 Mj), A 
transit time combined with the supposed period is given 
in Tabled 

Although this event is at a period approximately 4- 
times that of the known planet, it is not expected that 
this event is related to an alias of the data due to the 



TABLE 3 
Transit candidate parameters 



parameter value 



P 14.3053d 

R p 0.261Rj (2.92R e ) 

i 88.8° 

<j> 0.871 

Ephemeris Minimum (JD-2451545) 2051.4554 

Ax 2 %/A x 2 _% 0.893 

Mass (assuming p RJ 3000 kg m -3 ) 0.043 Mj (13.6 M s ) 



gaps associated with removing the transits of the known 
planet. The EEBLS algorithm is dissimilar to a Fourier 
transform in this regard, as gapped data of a certain 
period do not induce significant aliases at that period. 
Thus removing the transits of the known planet should 
have fully removed all signal at the period, harmonics, 
and subharmonics of the known planet. This has been 
confirmed by limited numerical tests using data with the 
same time-sampling as the data used in this applica- 
tion. A significant signal is not observed at harmonics, 
sub-harmonics of the known planet's orbital period sta- 
tistically often in these tests. This putative signal at 
P w 14.305g? is therefore intrinsic to the data. Evidence 
for a transit at this period is marginal, but additional 
MOST photometry of the HD 209458 system should con- 
firm or disprove its existence. 

The most significant brightening event was one ob- 
served with a period of approximately P « 8.96d It 
is likely statistical in nature, and thus is not expected to 
be related to any specific astrophysical process. 

As an additional sanity check this transit search 
method was applied to the current data-set without the 
removal of the transits of the known planet, HD 209458b. 
As expected, the routine correctly uncovers the transit 
of HD 209458b to a high deg ree of accuracy. The actual 
parameters of HD 209458b (jKnutson et all 12007m were 
recovered within 0.0004% and 0.4% in period and phase, 
respectively, while the inclination angle and radius were 
recovered within 1.5°, and 0.02i?j of the actual parame- 
ters. 

6. DISCUSSION 

MOSTs 2004 and 2005 observations of HD 209458 
have been searched for evidence of other exoplanets in 
the system. A transit-search routine has been adapted 
to search for more realistic limb-darkened transits to take 
advantage of the precise, near-continuous photometry re- 
turned by MOST and the fact that the stellar character- 
istics of the star, HD 209458, are well established. Monte 
Carlo statistics were generated using this routine and in- 
dicate that this routine in combination with the afore- 
mentioned MOST data has placed unparalleled limits on 
the size of transiting bodies that have been ruled out 
in this system for a range of periods, and inclination 
angles. In the most optimistic case of edge-on transits 
planets with radii greater than 0.20 Rj (~2.2 i?®) to 
0.36 Rj (^4.0 i?®) with periods from half-a-day to two 
weeks, respectively, have been ruled out with 95% con- 
fidence through this analysis. Specifically, we have been 
able to rule out transiting planets in this system with 
radii greater than those given in Figure [4] This work 
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has constrained theories that invoke smaller exoplanets 
in the system to explain the anomalously large size of 
HD 209458b, and those that predict "hot Earths" due to 
the inward migration of HD 209458b. 

6.1. Mass constraints on other Exoplanets in the system 

The above limits on the smallest planetary radii ex- 
cluded by MOST for transiting planets allow some es- 
timates on the type of planets that are being excluded. 
Planet radii change with time as the planet cools, so 
we make use of the fact that the age of the HD 209458 
system is determined to be about 5 Gyr with an uncer- 
tainty of 1.0 to 1.5 Gyr (jCodv fe Sasselovll2002ft . This 
means that we can safely refer to cooling models similar 
to the planets in the solar system. One exception would 
be the consideration of possible planets in orbits smaller 
than that of HD209458b, where tidal heating will lead to 
larger radii per given mass and composition. 

Our radius limit of 2.2-4.0 i?© already places us in 
the realm of super-Earth and hot Neptune planets. The 
range of bulk compositions could encompass iron-rich 
cores (super-Mercuries, ~ 80% Fe or Feo.s(-F 1 eS')o.2 
core), super-Earths (less than 50% Fe core), water- rich 
super-Earths (more than ~ 10% H 2 0), and water giants 
(hot Neptunes). Internal st ructure models f or such plan - 
ets have been published by I Valencia et al.1 (|2006l l2007f ). 
Hot Neptunes will have structures similar to Neptune 
(17.2 M e , 3.9 R®) and Uranus (14.4 M e , 4.0 R m ). Plan- 
ets with significant amount of water should be referred to 
as 'ice planets', because water ices VII & X would form 
at high pressures of a few GPa. At partial pressures ex- 
ceeding ~ lOOMPa ice planets are unlikel y to evaporate 
(in 5 Gyr) in orbits larger than 0.1 AU (Va lencia et al.l 
120071 ). while a pure H+He planet might (|Yellell2004h ~ 

In summary, our radius limit cannot exclude super- 
Mercuries of any Fe content, as their masses would have 
to exceed that of HD209458b itself. It can exclude some 
super-Earths, with masses larger than ~ 8 Mg. It ex- 
cludes water-rich super-Earths with masses larger than 
~ 6 M© for 50% water content, and ~ 5.5 M© for ice 
planets. The latter two types of planets will have radii 
that are very sensitive to stellar heating (hence their or- 
bit) and tidal heating. 

HD209458 is a system that has experienced migration 
and recent sim u lation s for such planetary systems by 
iRavmond et al.l (|2006h have confirmed previous expec- 
tations that water-rich Earths and super-Earths could 
often end up in orbits of 0.05 to 0.3 AU. Our transit 
search with MOST rules out super-Earths and hot Nep- 
tunes from 0.01 to 0.12 AU for HD209458 for edge-on in- 
clinations. Planets with inclination angles substantially 
different then edge-on (see Figured]) cannot be ruled out 
with this photometry, except for very short periods The 
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radius limits presented here, combin ed with the mass- 
limit s from transit-timing analyses (|Miller-Ricci et all 
[20061: lAgol fe Steffenl 12007). places very firm constraints 
on the size and masses of bodies that could still reside 
in close orbits to HD 209458. Reobservations of this sys- 
tem with MOST should be able to further reduce the size 
of planets that could remain undetected in this system, 
and further constrain theories invoking these planets in 
nearby orbits to HD 209458b. It will be very interesting 
to use MOST and obtain similar results for HD189733 
and HAT-P-1 and constrain such theoretical work. 
6.2. Radius and Mass constraints on Trojans 

We can also place a limit on the size of Trojans that 
may be leading or trailing HD 209458b, assuming the 
Trojan consistently transits HD 209458 with an orbital 
inclination close to that of the known pl anet (87°). We 
use th e definition of Trojans given by iFord fe Gaudi I 
(2006) of objects occupying the L4 or L5 stable La- 
grangian points. Trojans lagging or leading HD 209458b 
with a radius above R — 0.30 Rj (3.4 i?©) should have 
been detected with this analysis with 95% confidence. 
The 99% confidence contour is only slightly larger, R — 
0.31 Rj (3.5 i?©). If these putative planets were super- 
Earth analogues (p k, 3000 kg m~ 3 ) the corresponding 
mass limits would be M = 21M© and M = 23M© at 
95% and 99% confidence, respectively, while for hot Nep- 
tune analogues (p « 1600 kg m~ 3 ) the mass limits would 
be M = 11M© and M = 12M ffi . Depending on the 
structure of planet that is assumed, these mass limits are 
comparable or mo derately worse than th e 13 M© 99.9% 
confidence limit of lFord fe Gaudi] (p006l ) for HD 209458. 
However a detailed study using this MOST photometry 
of HD 209458 should be able to drastically improve this 
limit. It should be noted that lFord fe Gaudi I (|2006l ) cast 
doubt on the hypothesis that a hypothetical Trojan will 
consistently transit HD 209458. They note that if the 
hypothetical Trojan has a vertical libration amplitude 
greater than approximately 9° it will not consistently 
transit HD 209458. 
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Fig. 1.— The 2004 (A) and 2005 (B) HD 209458 MOST data-sets following the initial reduction and the automatic filtering step. The 
same data are shown binned in 30 minute intervals at a different vertical scale for 2004 (C) and 2005 (D). 
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Fig. 2. — EEBLS spectrum of the total MOST HD 20 9458 data-set, plotted in frequency (top), and period (bottom). The 4 candidates 
that passed the transit selection criteria as discussed in i|3.3l are shown by the solid vertical lines. In the frequency plot the lines marking 
the two candidates with the lowest frequencies appear as one. 
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Fig. 3. — A histogram showing the improvement of the transit over the anti-transit model (Ax 2 %/Ax 2 _%) for all Monte Carlo candidates 
that are spurious (when the transit recovered by the transit routine is not the inserted transit). The threshold level for Ax 2 %/Ax:l% 
above which a transit can be considered significant is determined via these statistics as it is set at a level, Ax 2 %/Ax?_% ~ 1.01, that would 
rule out 95% of these spurious signals. The solid-vertical line thus marks this threshold value. 
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Fig. 4. — Confidence limits in transit detection as a function of planet radius and planetary orbital period, for different orbital 
inclinations, based on Monte Carlo statistics. The crosses represent the radii and periods at which synthetic transits were inserted in the 
data. The dotted, thick-solid, and thin-solid lines represent the 99, 95 and 68% confidence contours, respectively. The thick dot-dashed 
curve represents the 68% confidence contour if the criterion A% 2 %/Ax?_% > 1.01 was not used. The near vertical dot-dash line in the 
later panels indicates the maximum period that produces a transit for that given inclination angle. Note the logarithmic period scaling on 
the x-axis. 110 phases were inserted for the 90° inclination angle case, while 65 phases were inserted for the other inclination angles. 
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Fig. 5. — Likelihood of spurious transit detections returned by the Monte Carlo statistical analysis. The dotted, thick-solid, and 
thin-solid lines represent the 10, 5 and 1% spurious signal contours, respectively. The format of the figure is otherwise identical to Figure 
[4] Note that the spurious signals occur in the intermediate regions of the period-radii space of interest, where it is not guaranteed that the 
correct transit will be recovered, but the inserted transit still causes significant deviations to the light-curve. 
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Fig. 7. — The best transit candidate as identified by our analysis of MOST's HD 209458 2004 and 2005 data-sets, but with marginal 
significance. The thick red-line represents the transit model, while the thin green-dashed line represents the constant brightness model. 
At top the data unbinned in 2004 (yellow) and the first (cyan), second (magneta), and third (blue) putative transit in 2005, and binned 
(black). Middle: the binned data only, at a different vertical scale. At bottom, an expanded portion of the phase diagram, with the binned 
data from the various transits given by the same colours as above. As this candidate failed the improvement in transit over anti-transit 
criterion, Ax 2 %/Ax?_% > 1.01, we report it as a possible, but unlikely candidate. The period and radius of this putative planet would 
be approximately 14.3d and 0.26 Rj (2.9R($). 



